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ABSTRACT 

Using ultra-violet absorption-lines, we analyze the systematic properties of the warm ionized phase 
of starburst-driven winds in a sample of 39 low-redshift objects that spans broad ranges in starburst 
and galaxy properties. Total column densities for the outflows are ~10 21 cm -2 . The outflow velocity 
( v out) correlates only weakly with the galaxy stellar mass (M*), or circular velocity (v C j r ), but strongly 
with both SFR and SFR/area. The normalized outflow velocity (v out /u cir ) correlates well with both 
SFR/area and SFR/M*. The estimated outflow rates of warm ionized gas (M) are ~ 1 to 4 times the 
SFR, and the ratio M/SFR does not correlate with v ou t- 
We show that a model of a population of clouds accelerated by the combined forces of gravity 
and the momentum flux from the starburst matches the data. We find a threshold value for the 
ratio of the momentum flux supplied by the starburst to the critical momentum flux needed for the 
wind to overcome gravity acting on the clouds ( R cr it )• For R cr it > 10 (strong-outflows) the outflows 
momentum flux is similar to the total momentum flux from the starburst and the outflow velocity 
exceeds the galaxy escape velocity. Neither is the case for the weak-outflows ( R cr it < 10). For the 
weak-outflows, the data severely disagree with many prescriptions in numerical simulations or semi- 
analytic models of galaxy evolution. The agreement is better for the strong-outflows, and we advocate 
the use of R cr it to guide future prescriptions. 

Subject headings: galaxies: starbursts — galaxies: ISM —galaxies: intergalactic medium—galaxies: 
evolution—galaxies: kinematics and dynamics 


1. INTRODUCTION 

In the standard picture, a galaxy will form when gas 
is able to reach high enough densities to cool, sink to 
the center of its dark matter halo, and form stars. What 
happens next depends on the interplay between a host of 
complex physical processes. In particular, massive stars 
produce a copious supply of ionizing radiation, and then 
quickly end their lives as supernovae. These processes 
ionize, heat, and chemically-enrich their surroundings. 
Unfortunately, a robust theoretical understanding of stel¬ 
lar feedback has proved elusive. The physical processes 
that control star formation and the resulting feedback 
operate on scales that are still well below the resolu¬ 
tion of any cosmological simulation (‘sub-grid physics’). 
These feedback processes are normally parametrized in 
both numerical simulations and semi-analytic models us¬ 
ing prescriptions based on simple theoretical arguments 
and/or empirical relations (see Somerville & Dave 2015 


for a recent review). 

Galactic winds flowing from intensely star-forming 
galaxies are an especially important form of feedback. 
They are driven by the energy and/or momentum sup¬ 
plied by massive stars and their evolutionary descendants 
(e.g. Chevalier & Clegg 1985; Murray et al. 2005). 
Winds like these are expected to alter the structure and 
evolution of the galaxy, possibly slowing down the rate 
at which it can subsequently form stars. They can also 
expel metal-enriched material into the circum-galactic 
medium and even into the inter-galactic medium. These 
winds are complex multi-phase phenomena, with out¬ 
flowing material spanning the range from cold molecu¬ 
lar gas revealed in the mm and sub-nun regime to very 
hot gas studied in X-rays and the relativistic magnetized 
material seen in synchrotron emission in the radio (see 
Veilleux et al. 2005 and references therein). Thus, the 
full range of the wind properties and effects can only be 
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studied in relatively nearby galaxies. 

In the contemporary universe, strong winds are only 
observed in galaxies undergoing intense bursts of star 
formation. In contrast, at high redshift, winds are es¬ 
sentially ubiquitous in typical star-forming galaxies (e.g. 
Shapley 2011). Ideally then, the insights provided by de¬ 
tailed observations of winds in nearby starbursts should 
be combined with the demographic information about 
winds at high redshift, in order to understand the role of 
galactic winds in the evolution of galaxies and the inter- 
galactic medium. 

The great majority of information about galactic winds 
at high redshift has come from observations of interstel¬ 
lar absorption lines in rest-frame ultraviolet (UV). These 
lines trace outflowing warm neutral and ionized gas (e.g. 
Shapley 2011; Steidel et al. 2010). 1 They can also 
be observed in local starburst galaxies using space-based 
UV observations. The goal of the present paper is to 
use UV spectroscopic observations of a well-chosen and 
well-characterized sample of low-z starburst galaxies to 
determine which properties of the starburst and the sur¬ 
rounding galaxy correlate best with the principal prop¬ 
erties of the outflows. Our hope is that these empirical 
correlations will not only lead to a better understand¬ 
ing of galactic winds, but can also be used to inform the 
sub-grid treatment of galactic winds in future numerical 
simulations. 

In section 2 we describe the two samples of starbursts 
and summarize the analysis of the UV spectra. In sec¬ 
tion 3 we present empirical results on the correlations 
of measured outflow velocity and estimated mass and 
momentum outflow rates with the starburst and galaxy 
properties and compare our results with previous stud¬ 
ies. In section 4 we compare our observations with 
simple analytic models for momentum-driven winds and 
discuss the implications for models and simulations of 
galaxy evolution. We then summarize our results in sec¬ 
tion 5. Throughout the paper we adopt a cosmological 
model with Hq = 70 km s -1 Mpc -1 , Ha = 0.7, and 
^Matter 0.3. 

2. SAMPLE SELECTION & DATA ANALYSIS 
2.1. The Samples 

We use two samples of low-redshift (z < 0.2) starburst 
galaxies in the analysis presented in this paper. The first 
is a sample of 19 galaxies observed with the Far Ultravi¬ 
olet Spectroscopic Explore (FUSE), and was previously 
investigated by Grimes et al. (2009: hereafter G09) and 
Heckman et al. (2001). The second is a sample of 21 Ly¬ 
man Break Analogs (LBAs) observed with the Cosmic 
Origins Spectrograph (COS) on the Hubble Space Tele¬ 
scope (HST), and was previously investigated by Alexan- 
droff et al. (2015: hereafter A15). These papers contain 
descriptions of both the sample selection and the data 
analysis. We refer the reader to these papers for the 
details, but as a matter of convenience we give a brief 
summary here. 

1 In this paper we use the term ‘wind’ to describe the full 
multi-phase phenomenon. We use the term ‘hot wind fluid’ to 
describe the outflowing thermalized ejecta from supernovae and 
stellar winds. A combination of this hot wind fluid and stellar ra¬ 
diation interacts with the interstellar or circum-galactic medium to 
produce the outflowing gas seen in absorption (which we will term 
the ‘outflow’). 


The COS LBA sample was drawn from the large par¬ 
ent sample derived from a cross-match of GALEX GR6 
and SDSS DR7 (Overzier et al. 2015). The criteria 
used to define the parent sample were a K-corrected 
GALEX-based FUV luminosity Lpuv > 10 103 Lq and a 
GALEX based FUV effective surface-brightness Ifuv(= 
l/ 2 LFuv/ 7rr 5 o) > 10 9 Lgkpc -2 , where r 50 was the 
seeing-deconvolved SDSS u-band half-light radius (Heck¬ 
man et al. 2005; Hoopes et al. 2007). These criteria were 
designed to select objects with luminosities and surface- 
brightnesses similar to typical Lyman Break Galaxies 
(LBGs) at z ~ 3 to 4 (e.g. Shapley 2011). Subsequent 
observations have shown that the LBAs closely resem¬ 
ble LBGs in all salient properties: stellar mass, SFR, 
metallicity, galaxy size and morphology, gas kinematics, 
and gas and dust content (Hoopes et al. 2007; Basu- 
Zych et al. 2007; Overzier et al. 2010,2011; Gongalves et 
al. 2010,2014). The specific targets observed with COS 
were selected based on maximizing the estimated FUV 
flux within the COS aperture. These data include COS 
G130M and G160M spectra for each galaxy, and were 
obtained as part of the HST Cycle 17 Program 11727 
and Cycle 20 Program 13017 (PI: T. Heckman). 

The FUSE starburst sample was drawn from the 
archival collection of data on low-redshift star-forming 
galaxies. The specific members of this sample were se¬ 
lected for analysis solely on the basis of data quality (S/N 
> 4 in 0.078 A bins in the galaxy continuum in the LiF 
1A FUSE band). These galaxies are broadly representa¬ 
tive of the local population of UV-bright starburst galax¬ 
ies. Three members of the FUSE sample (VV 114, Haro 
11, and Mrk 54) meet the LBA criteria (Grimes et al. 
2006,2007; Hoopes et al. 2007). 

We have selected these two samples for several reasons. 
First, the spectroscopic aperture used is large enough to 
encompass the majority of the far-ultraviolet light emit¬ 
ted by the starburst. More quantitatively, the size of the 
COS aperture exceeds the UV half-light diameters (see 
below) of the galaxies by an average factor of 8 and by a 
minimum of factor of 3. The size of the FUSE aperture 
exceeds the UV half-light diameters (see below) by an 
average factor of 4.4. Only in the case of M 83 is the 
FUSE aperture smaller than the galaxy UV half-light di¬ 
ameter. This insures that we are measuring the global 
properties of the starburst in nearly every case. 

Secondly, these two data sets are complementary. As 
noted above, the COS LBA sample is representative 
of typical star-forming galaxies at high redshift. The 
FUSE starburst sample extends to galaxies of lower stel¬ 
lar mass, star-formation rates, and metallicity. This is 
important, since the goal in this paper is to quantify the 
dependence of the wind properties on these basic param¬ 
eters. Third, both the FUSE and COS data have point 
source spectral resolutions (R ~ 21,000 and 18,000 re¬ 
spectively). Even for extended targets, the effective res¬ 
olution is sufficient to resolve the interstellar absorption¬ 
lines in these galaxies (unlike IUE, HUT, or HST/STIS). 
It is also important to note that we have undertaken the 
analysis of both sets of data and have been uniform in the 
way we have measured the properties of the interstellar 
absorption-lines used to probe galactic winds. 

This sample complements previous similar investiga¬ 
tions. With the exception of the optically-selected sam- 
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pie in Chen et al. (2010), the previous analyses of compa¬ 
rably large samples of low-redshift galaxies have focused 
on dusty far-IR-selected galaxies (e.g. Heckman et al. 
2000; Rupke et al. 2005; Martin 2005; Sato et al. 2009). 
All these low-z studies have used the optical Na I D dou¬ 
blet to probe the outflows. Our sample is essentially 
UV-selected, making it a better match to typical high-z 
star-forming galaxies (e.g. LBGs), and giving us access 
to diagnostic interstellar lines in the far-UV spanning a 
range in ionization state and optical depth. 

There have also been a number of investigations of in¬ 
termediate redshift (z ~ 0.5 to 1.5) star-forming galaxies 
(e.g. Weiner et al. 2009; Erb et al. 2012; Kornei et al. 
2012; Martin 2012; Rubin et al. 2014; Bordoloi et al. 
2014). These have typically used stacked spectra of the 
Mg II 2800 doublet and/or several near-UV Fe II lines. 
Our sample spans a much larger range in star-formation 
rate and galaxy mass, giving us better leverage on the 
influence these properties have on those of the outflows. 
We also have adequate signal-to-noise and spectral reso¬ 
lution to be able to characterize the outflows in individual 
galaxies. 

The principal properties of our merged COS and FUSE 
sample are listed in Table 1. The starburst half-light radii 
(r*) for the LB As were measured using HST UV images 
(A15), and from the available UV images described in 
G09 for the FUSE sample. The stellar masses (M*) for 
the LBAs were taken from the MPA-JHU galaxy catalog 
for SDSS DR7. 2 The stellar masses for the FUSE sam¬ 
ple were based on 2-MASS K-band luminosities and the 
K-band mass/light ratio for strongly star-forming galax¬ 
ies derived by Bell & De Jong (2001). Star-formation 
rates (SFR) were derived using a combination of GALEX 
FUV, IRAS and Herschel far-IR, and WISE and Spitzer 
mid-IR photometry following the prescriptions in Ken- 
nicutt & Evans (2012), as described further in A15. In 
all cases, a Kroupa/Chabrier Initial Mass Function was 
adopted for masses and SFRs. We have then defined an 
effective star-formation rate per unit area as SFR/2nrl. 
The metallicities (O/H) are based on the 03N2 strong¬ 
line estimator derived by Pettini & Pagel (2004), using 
SDSS spectra for the LBA sample, and literature spectra 
for the FUSE sample. 

We have determined a characteristic internal veloc¬ 
ity parameter for each galaxy. Gongalves et al. (2010) 
showed that both large-scale rotation and turbulent mo¬ 
tions were significant in the kinematics of the interstellar 
medium of LBAs, and this is also seen in more typical 
local starburst galaxies (e.g. Lehnert & Heckman 1996). 
We therefore use the kinematic parameter S (Weiner et 
al. 2006; Kasim et al. 2007) to characterize each galaxy: 
S 2 = 0.5v 2 ot + <j 2 as (where v ro t and a gas represent the 
observed rotation velocity and velocity dispersion of the 
ISM). We adopt a value based on the good fit between S 
and M* for low-redshift star-forming galaxies shown in 
Simons et al. (2015): log S = 0.29 log M* — 0.93, where 
S is in km sec -1 , and M* is in solar masses. We find 
that this relationship is a good fit to the data for the 16 
LBAs in Gongalves et al. (of which 6 are in the COS 
sample): the mean residual of the (observed minus fit) 
value for log S is -0.01 dex, with a standard deviation of 
0.10 dex. In the rest of the paper we use the value of S 

2 Available from http://www.mpa-garching.mpg.de/SDSS/DR7/ 


from this fit to estimate an effective circular velocity for 
each galaxy: v C j r = y/2S (where the y/2 factor insures 
that v C i r = v ro t in the limit v ro t » &gas)- 3 

2.2. Analysis of the UV Spectra 

We refer the reader to G09 and A15 for detailed dis¬ 
cussions of the data reduction for the FUSE and COS 
data respectively. In this paper we use the reduced spec¬ 
tra reported in these papers to measure two quantities: 
the column densities and velocity of the outflow based 
on UV interstellar absorption-lines. 

2.2.1. Column Densities 

The strong interstellar lines in these spectra are sat¬ 
urated (highly optically-thick) and therefore cannot be 
used to derive reliable ionic column densities. Instead, 
we will use optically-thin (and therefore weak) lines. In 
order to measure these weak features in the COS sample, 
we have created a high signal-to-noise stacked spectrum. 
The procedure is discussed in A15 and the resulting stack 
is shown that paper. From this stack we have used the 
following weak lines to trace both the neutral and ionized 
gas: S III 1012, O 1 1039; S IV 1063, Fe II 1145, 1608; Al 
II 1670. The FUSE data are higher in signal-to-noise, so 
that the following weak lines were measured by G09 for 
individual galaxies: S III 1012, Si II 1021, OI 1039, S IV 
1063, and Fe II 1097, 1145. 

Ionic column densities were obtained by integrating the 
apparent optical depths of these lines over the velocity 
range covered by the absorption (Sembach et al. 2003; 
G09), using oscillator strengths from the on-line NIST 
Atomic Spectra Database. 4 These column densities can 
be used to estimate total column densities for the neutral 
and warm ionized gas. The neutral phase is traced most 
naturally by O I. However, S II, Si II, Fe II, and Al 
II are the dominant ionic species of these elements in 
the warm neutral atomic (HI) phase of the ISM in the 
starbursts since they all have ionization potentials less 
than a Rydberg. We use the sum of the S III and S IV 
column densities to trace the ionized gas. 

For the COS stacked spectrum we 
measured the following column densities 
(cm -2 ) : Noi < 10 15 - 7 , N FeII = 10 146 , N A m = 10 13 - 4 . 
We have converted these into HI column densities using 
the mean O/H abundance derived from the nebular 
emission-lines for the sample (0.5 solar) and assuming 
solar O/Fe and O/Al abundance ratios. These yield 
values for N HI of < 10 193 , 10 194 , and 10 19 3 cm -2 
based on the O I, Fe II, and Al II lines respectively. The 
column density derived from O I is most reliable (least 
sensitive to dust depiction or ionization corrections). 
The estimated column densities are consistent with 
the lack of damping wings on the Lya absorption-lines 
in the COS data for our LBAs. The same exercise 
yields column densities of Nsm = 10 15 3 cm -2 and 

3 We also estimated values for the halo virial velocity (v v i r ) 
using the scaling relation between M* and halo virial mass from 
Behroozi et al. (2010). Over the range in stellar mass where their 
relationship is calibrated ( log M* > 9.0), the values for v v i r agree 
with our estimates of v c i r to within 0.1 dex. We use v c i r in this pa¬ 
per because its calibration has been more directly determined from 
data for strongly star-forming galaxies, and because it is calibrated 
over the full range of M* we are probing. 

4 http: //www. nist. gov/pml/data/asd. cfm 
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Fig. 1.— We compare the absorption-line profiles of five strong 
interstellar transitions, extracted from the stacked, smoothed 75 
km s —1 FWHM) high signal-to-noise spectrum of the COS sam¬ 
ple described in Alexandroff et al. (2015). The continuum was 
normalized to unit flux locally for each of the lines in the figure. 
The figure shows that the profile shape of the Si III 1206 line that 
we use to measure outflow velocities in the COS sample (shown in 
purple) is very representative of the profile of the N II 1084 line (in 
blue) that we use for the FUSE sample, and to the profiles of the 
other strong lines (Si II 1260 in yellow, C II 1334 in green, and 
Si IV 1393 in orange). The velocities on the x-axis are measured 
relative to the galaxy systemic velocity. 

Nsiv = 10 15 - 2 cm -2 . Sulfur is not significantly depleted 
onto dust grains, and we do not expect a significant 
amount of S V. The combined column density of these 
two ions corresponds to Nhii = 10 20 7 cm -2 . 

For the FUSE data we can use the same approach but 
derive column densities for individual galaxies, convert¬ 
ing ionic to Hydrogen column densities using the metal- 
licities for each galaxy listed in Table 1. We see no ev¬ 
idence for significant systematic variations in H column 
densities from galaxy to galaxy, so we simply calculated 
the mean values for the FUSE sample: N hi = 10 19 ' 6 and 
N hi i = lO 20 6 cm -2 . These are similar to the values for 
the COS sample. 

2.2.2. Outflow Velocities 

It is usually thought to be straightforward to measure 
outflow velocities using the observed absorption-line pro¬ 
files. However, as emphasized by Prochaska et al. (2011) 
and Scarlata & Panagia (2015), the UV absorption-lines 
are resonance lines and as such, each upward transition 
resulting in absorption should be followed by a downward 
radiative transition resulting in emission. This emission 
can infill the absorption-line, and can in principle have 
a significant effect on the net observed profile. Both of 
the above papers show that the fine-structure transitions 
that are associated with many of the commonly-used 
resonance lines can be used to deduce the properties of 
the infilling resonance line emission. These effects have 
been investigated observationally in several recent papers 
(Martin et al. 2012; Erb et al. 2012; Kornei et al. 2012; 
Rubin et al. 2014). 

A similar approach to these observational papers has 
been adopted by A15, who have explicitly shown that the 
effect of infilling has a negligible effect on the properties 
of the outflows inferred from the observed absorption¬ 
line profiles in the COS sample. We have inspected the 
individual spectra of the FUSE sample, looking for the N 


II* 1085.7 and N III* 991.5 fine structure emission-lines 
(associated with the N II 1084.0 and N III 989.8 reso¬ 
nance lines). No significant detections were found, and 
therefore infilling will not be significant in these galaxies 
either. 

Given that the ionized gas is the dominant phase in 
the outflows (see above), we have chosen to character¬ 
ize the outflow velocity using the strongest relevant lines 
available. For the COS data, this is the Si III 1206 line 
(A15). In Figure 1 we show line profiles for several strong 
interstellar absorption lines taken from the high signal- 
to-noise, smoothed (R ~ 4000) stacked spectrum of the 
COS sample described in A15 (N II 1084, Si III 1206, Si 

II 1260, C II 1334, Si IV 1393). It is clear that the Si 

III 1206 line is representative of the profiles of the other 
strong lines. For the FUSE data, the strongest line trac¬ 
ing ionized gas is CIII 977, but the line measured with 
the highest S/N ratio is Nil 1084 (G09). We take the 
un-weighted mean of the measured outflow speeds using 
these two lines. As seen in Figure 1, the N II absorption¬ 
line in the stacked COS spectrum has the same properties 
as the Si III line, so there should be no systematic dif¬ 
ferences between the outflow velocities measured in the 
FUSE vs. COS samples. 

In all cases we use a non-parametric measure of the out¬ 
flow velocity based on the flux-weighted line centroid de¬ 
fined relative to the systemic velocity of the galaxy (G09; 
A15). We emphasize that, absent a good model for the 
radial dependence of the mass-weighted distribution of 
the outflow velocity field, any specific definition of a char¬ 
acteristic outflow velocity is somewhat arbitrary. While 
a different outflow velocity from ours could be adopted 
(for example, the maximum value where the line meets 
continuum, or the velocity enclosing the 80% red-most 
portion of the integrated line profile), such parameters 
are usually more difficult to measure (more subject to 
uncertainties in the continuum placement and more de¬ 
pendent on the signal-to-noise in the data). Our choice 
is based on ease and robustness of measurement, and 
does not imply some specific model for the outflow. In 
any case, using a different definition would not affect the 
qualitative sense of our results. 

The line centroids are significantly blue-shifted in all 
21 of the LBAs with COS data and in 15 of the 18 FUSE 
starbursts. This underscores the ubiquity of outflows in 
local starbursts. These outflow velocities are given in 
Table 2. 

3. EMPIRICAL CORRELATIONS 
3.1. Outflow Velocity 

We begin by examining the correlation of the outflow 
velocity (v out ) with the principal properties of the galaxy 
and its starburst. Throughout the remainder of the pa¬ 
per we will use the Kendall r test to assess the statistical 
significance of each correlation. As shown in Figure 2, 
we see only a weak correlation with M* (and hence with 
v C i r ). There is a statistically significant positive corre¬ 
lation with SFR, SFR/M* , and SFR/area. The last 
correlation is the strongest, and the impression is that 
the outflow velocity climbs rapidly with SFR/area but 
then saturates at about 400 km s -1 for SFR/area > 6 
M 0 yr -1 kpc -2 . We will return to this result in sec¬ 
tion 4 below where we compare the data to predictions 
of idealized wind models. 
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Fig. 2. — The log of the outflow velocity is plotted as a function of the basic properties of the starburst galaxies. Clockwise from the 
upper left (a,b,c,d). The upper left panel (a) shows that there is only a weak correlation between the outflow velocity and the galaxy 
circular velocity. The diagonal lines show v ou t = 10 v c i r , v ou t = v cvr . v ou t = 0.1v c i r . The ratio v ou t/v c i r spans a range of over two 
orders-of-magnitude. The label on the upper axis shows the corresponding values of the galaxy stellar mass (see text). The upper right 
panel (b) shows a significant correlation between SFR, and v ou t, as has also been seen in other previous studies (see text). The bottom two 
panels show the correlation with two forms of normalized SFR: SFR/area (c) and SFR/M* (d). Both correlations are significant, but the 
correlation with SFR/area is stronger. The crosses represent the typical uncertainties (see Table notes). The blue and green points show 
the strong- and weak-outflows. See section 4.1 and Fig. 9a. 


The correlation between outflow velocity and star- 
formation rate (albeit with significant scatter) has been 
previously seen in samples of low-redshift starbursts 
(Heckman et al. 2000; Martin 2005; Rupke et al. 2005; 
Martin et al. 2012), with a form that is qualitatively 
consistent with our Figure 2b. Studies at higher-redshift 
have been more mixed, but this most likely reflects the 
rather small range in star-formation rate that has been 
probed (Weiner et al. 2009; Erb et al. 2012; Kornei et al. 
2012; Rubin et al. 2014; Bordoloi et al. 2014). There has 
also been disagreement about whether the outflow veloc¬ 
ity correlates with SFR/area (Chen et al. 2010; Kornei 
et al. 2012; Rubin et al. 2014), but again these studies 
have covered relatively small ranges in this parameter 
compared to our sample. 

The possible dependence of outflow velocity on either 
M* or v c i r has been examined by many of these stud¬ 
ies. Most have found either no correlation, or only a 


weak one (Heckman et al. 2000; Rupke et al. 2005; 
Erb et al. 2012; Martin et al. 2012; Rubin et al. 2014; 
Bordoloi et al. 2014). The strongest correlation was 
found by Martin (2005). Ours is the only investigation 
in which the dynamic range in these quantities was sim¬ 
ilar to hers. Unfortunately, her sample only consisted 
of nine galaxies. The galaxies with low stellar masses 
and circular velocities in her sample did not contain the 
population of objects with relatively high outflow veloc¬ 
ities ( v ou t ~ 10 2 - 5 km sec -1 ) and low circular velocities 
( v C ir < 10 2 km sec -1 ) seen in our Figure 2b. We empha¬ 
size that we only poorly sample the population of galax¬ 
ies with M* < 10 8 Mg ( v C i r < Sbkms -1 ) and cannot say 
whether the trends we see extend into this range. 

We can also consider the normalized outflow velocity 
(v out/v C ir) ■ As seen in Figure 3, this ratio spans over 
two-orders-of-magnitude (~< 0.1 to 10). There are sig¬ 
nificant positive correlations between the ratio and nor- 
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plotted as a function of the specific SFR (SFR/M»). Both plots show strong correlations. The crosses represent the typical uncertainties 
(see Table notes). The blue and green points show the strong- and weak-outflows. See section 4.1 and Fig. 9a. 
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log N h 

Fig. 4.— A histogram of the total Hydrogen column density for 
the merged COS plus FUSE sample, estimated using the amount 
of the effective UV dust extinction, the Calzetti et al. (2000) star- 
burst attenuation law, and an assumption of a fixed (Milky Way) 
ratio of dust to metals in the ISM. The typical column density 
of 10 21 cm' 2 agrees with estimates based on the optically-thin 
interstellar absorption lines in these galaxies. See text for details. 

malized measures of the star formation rate (SFR/area 
and SFR/M*). The correlation with SFR/area also 
shows a suggestion of a possible saturation in normalized 
outflow velocity at v ou t ~ 2 to 4 v C i r above SFR/area 
~ 10 M 0 yr -1 kpc~ 2 . It is important to emphasize 
that outflow velocities this high are well in excess of the 
galaxy escape velocities. Moreover, these velocities refer 
to the centroid of the absorption-line, while outflowing 
gas is observed at velocities at least twice this high (e.g. 
Figure 1). 

3.2. Outflow Rates 

In order to turn the outflow velocities into outflow 
rates, several additional pieces of information are needed. 
More specifically, the mass outflow rate is given by 


M = ClN out < m > v out r out . Here, H is the solid angle 
occupied by the outflow, N oll t is the total column density 
of the outflow along the line-of-sight, < m > is the mean 
mass per particle, and r out is the characteristic radius 
of the location of the absorbing material in the outflow. 
Given the ubiquity of outflows in our sample, we will set 
Q = 47t (while acknowledging that our UV-bright sam¬ 
ple may be selectively missing galaxies viewed near the 
equatorial plane of a bi-polar outflow see Heckman et 
al. 2000; Chen et al. 2010; Kornei et al. 2012; Rubin 
et al. 2014). In the sections below we will describe the 
estimates of the column densities and radii of the winds. 

3.2.1. Estimated Column Densities & Wind Radii 

In section 2.2.1 we used weak optically-thin absorption¬ 
lines to estimate column densities in the outflows. To 
check whether the total inferred column densities are 
plausible, we have used an indirect method that is in¬ 
dependent of these estimates. We first defined an ef¬ 
fective far-UV opacity due to dust grains: Tdust = 
— ln(Luv,obs/ Luv,int )■, where the UV luminosities re¬ 
fer to the observed (transmitted) and intrinsic (inci¬ 
dent) quantities. We define huv,int = L UVi0bs + L IRito t, 
where L mjot is the luminosity between wavelengths of 
8 and 1000 /im and l-iuv,obs = 1-66 Lpuv is the luminos¬ 
ity shortward of 3000 A (Meurer, Heckman, & Calzetti 
(1999). See A15 and G09 for a description of the relevant 
data. The resulting values for Tdust are listed in Table 
1. The starbursts have a mean value Tdust = 1.5 with 
a dispersion of 1.0 (these objects are translucent in the 
UV). 

The Calzetti starburst dust attenuation law (Calzetti 
et al. 2000) then yields the column density of dust cor¬ 
responding to a given value of Tdust , and we convert this 
into a total Hydrogen column density assuming that all 
the galaxies have the same fraction of metals residing 
in dust grains (0.5) as in the Milky Way (see Matts- 
son et al. 2014 and references therein). The implied 
column densities are listed in Table 2 and shown in his¬ 
togram form in Figure 4. The distribution has a mean of 
N h = 10 210 cm -2 with a standard deviation of only 0.25 
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Log (SFR [Mg/yr]) 

Fig. 5. — The estimated outflow rate of ionized gas plotted as a 
function of the SFR (both in Mq year -1 ). The two quantities are 
well-correlated. The diagonal line shows a ‘mass-loading factor’ 
(M 0 ut/ SFR) of two, close to the median value for the sample. The 
cross represents the typical uncertainties (see Table notes). The 
blue and green points show the strong- and weak-outflows. See 
section 4.1 and Fig. 9a. 

dex. While this method is clearly a crude one, it is re¬ 
assuring that this column density is similar to the values 
derived more directly using the weak absorption lines. It 
also suggests that there is not a large range in total col¬ 
umn density within our UV-selected sample of galaxies. 
Higher column densities (few xlO 21 to 10 22 cm -2 ) have 
been inferred in the dusty outflows seen in Na I optical 
absorption-lines in far-IR-selected samples (Heckman et 
al. 2000; Rupke et al. 2005; Martin 2005). 

The radius to be used to calculate the outflow rate 
depends upon the assumed model for the absorbing gas. 
In the idealized case of a thin shell-like structure, it is 
simply the radius of that shell. For a mass-conserving 
outflow with v out independent of radius, the appropriate 
value to use is that of the starting point of the wind. 
For the classic case of an outflow of gas clouds driven 
by the ram pressure of a hot wind fluid created from 
thermalized supernova ejecta (Chevalier & Clegg 1985), 
this would be the sonic radius of the wind (roughly, the 
radius of the starburst). Since the correct value of the 
radius is uncertain, we will simply assume that it is some 
scale factor /3 times the starburst radius r*. 

3.2.2. Resulting Correlations 

Given the above, it is clear that the calculated outflow 
rates are uncertain, and should only be taken as roughly 
indicative estimates that are prone to systematic uncer¬ 
tainties. With that caveat, we list the estimated outflow 
rates in Table 2 and show the relationship between M out 
and the SFR in Figure 5. There is a significant corre¬ 
lation. The plot shown assumes a total wind column H 
density of 10 20 85 cm -2 (the mean of the estimates given 
above) and /? = 2. For these choices, the outflow rates 
typically range between 1 and 4 times the SFR (see Ta¬ 
ble 2). Similar results have been found previously for 
dusty low-z starbursts (e.g. Heckman et al. 2000; Rupke 
et al. 2005; Martin 2005), for intermediate redshift UV- 
selected galaxies (Bordoloi et al. 2014), and for higher- 
redshift Lyman Break Galaxies (Pettini et al. 2002). 


The ratio of M out /SFR is referred to as the mass¬ 
loading factor. In Figure 6 we show that the mass¬ 
loading factor shows no significant dependence on v ou t or 
r*, and only weak (inverse) dependences on M*, or v C i r . 
The apparent inverse correlation between M out /SFR 
with SFR must be treated with caution, since this a plot 
of y/x vs. x (which can induce a non-physical correla¬ 
tion). We will discuss the significance of these results in 
section 4 below. 

We can also estimate the momentum flux carried by 
the outflow: p out = Mv ou t- These values are listed in 
Table 2. We can now compare these values to the total 
momentum flux supplied by the starburst. This would be 
due to a combination of a hot wind fluid driven by the 
thermalized ejecta of massive stars (Chevalier & Clegg 
1985) and radiation pressure (Murray et al. 2005). Based 
on Starburst99 (Leitherer et al. 1999), for a typical star- 
burst population, the ratio of the momentum flux from 
these two respective sources is about 2.5. Summing them 
together yields a total momentum flux p* = 4.8 x 10 33 
SFR dynes, where the SFR is in M 0 year -1 . The re¬ 
sulting values are given in Table 2. In Figure 7 we find 
a strong but non-linear correlation between p out and p* 
(see also Rupke et al. 2005). We will discuss the signifi¬ 
cance of this below. 

4. DISCUSSION 

4.1. Comparisons to Simple Analytic Wind Models 

In this section we will compare the observed outflow 
velocities with predictions of several idealized analytic 
models for the outflow. 

4.1.1. Momentum-Driven Population of Clouds 

One simple model (e.g. Chevalier & Clegg 1985) is 
that the absorption-lines are produced by a population 
of clouds/filaments that are driven outward by the mo¬ 
mentum supplied by the starburst, as quantified above. 
In this picture, each cloud would be smaller than the 
starburst, and the observed absorption-line profile would 
reflect the distribution in radial velocity of the cloud 
population located along the lines-of-sight to the star- 
burst. Such a picture is consistent with optical imaging 
of the ionized emission-line gas seen in nearby starburst- 
driven outflows (e.g. Mutchler et al. 2007) and with 
high-resolution numerical simulations of such flows (e.g. 
Cooper et al. 2008, 2009). A model of this kind would 
require enough clouds along the lines-of-sight to the star- 
burst to produce the smooth absorption-line profiles that 
are observed in our galaxies (A15, G09). 

The outward force on a cloud with a cross-sectional 
area A c located a distance r from the starburst is: 

F out = A c p^/4Trr 2 (1) 

In addition to the outward force, a cloud with mass M c 
will also experience an inward force due to gravity: 

F in = M c u 2 ir r -1 (2) 

where v C i r = GM{< r)/r 1 ^ 2 and M{< r ) is the total 
mass interior to r. Writing M c as A C N C < m > (where 
N c is the cloud Hydrogen column density and < m > 
is the mean mass per H particle), we can then define a 
critical momentum flux such that the net force on a cloud 
located at r = r* is outward: 

Pcrit,c = 4-7T r*V c < m > v\ ir (3) 
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Fig. 6.— The mass-loading factor (M ou t / SFR) is plotted as function of the basic properties of the galaxy and its starburst. Clockwise 
from upper left (a,b,c,d). The two left panels (a,d) show weak inverse correlations of mass-loading with (top) SFR and (bottom) galaxy 
circular velocity and stellar mass. The two right panels (b,c) show that there is no correlation with either the outflow velocity (top) or 
the starburst radius (bottom). See text for details. The crosses represent the typical uncertainties (see Table notes). The blue and green 
points show the strong- and weak-outflows. See section 4.1 and Fig. 9a. 


In convenient units, p cr it,c is 10 33 9 dynes for N c = 
10 21 cm -2 , r* = 1 kpc, and v C i r = 100 km sec -1 . The re¬ 
sulting values for p cr it,c are listed in Table 2. In Figure 8 
we plot this critical momentum flux vs. the momentum 
flux available from the starburst. We then see that in 
95% of the cases (37/39), the starburst momentum flux 
exceeds the critical value for an outflow (often by over 
an order-of-magnitude). 

Taking this a step further, Figure 9 plots the normal¬ 
ized outflow velocity v 0 v.t/v C ir vs. p*/p cr it,c (hereafter 
Rcrit,c )• This shows a strong positive correlation. A 
closer examination suggests that v ou t/v C i r at first in¬ 
creases with increasing R cr it,c , but then saturates at 
about Vout/vtir ~ 3 to 5 for ratios of R cr it,c > 10. 
Crudely speaking, this would suggest three regimes: a 
‘no-outflow’ regime where p* < p C rit,a a ‘weak-outflow’ 
regime where p* ~ 1 — 10 p C rit,c, and a ‘strong-outflow’ 
regime, where p* > 10 p C rit,c- In this regime the out¬ 


flow exceeds the galaxy escape velocity. To illustrate 
their different properties, we have color-coded the strong- 
outflows (blue) and weak/no-outflows (green) in Figures 
2, 3, 5, 6, and 7. We will return to the implications of 
some of these differences in section 4.3 below. 

We now compare the predicted and observed outflow 
velocities. To do this, we adopt a simple model in which 
a given cloud starts at some initial radius at rest (we take 
r*) and is then accelerated by the combined force due to 
gravity and the starburst momentum flux. This yields 
the following expression for the velocity of the cloud as 
a function of its distance from the nucleus (/3r*): 

v c {p)/v cir = \/2[(l - ^ -1 )(i? c „ t , c ) - M/3)] 1/2 (4) 

The key importance of R cr it,c is clear in this expression. 
In the limit i? cr , itc >> 1, gravity becomes negligible and 
for (3 » 1 the expression for v c simplifies to that given 
in Chevalier & Clegg (1985). One possible explanation 
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Fig. 7. — The estimated observed rate of momentum transport 
in the outflow is plotted as a function of the rate of momentum 
supplied by the starburst (in dynes). The diagonal lines show the 
ratios p ou t/p* = 10, 1 , and 0.1. The forceful outflows ( p ou t > 10 34 
dynes) roughly carry the full amount of the momentum supplied 
by the starburst. The less forceful outflows carry typically only 
~ 10% of the available momentum. The cross represents the typical 
uncertainties (see Table notes). The blue and green points show 
the strong- and weak-out flows. See section 4.1 and Fig. 9a. 

for the observed spread of radial velocities of the clouds 
(as seen in the absorption-line profiles) wonld be a range 
in cloud column densities (since p C rit,c -V c )- 
Because the momentum flux acting on a cloud drops 
like 1 /r 2 , while the gravitational force in our assumed 
isothermal potential drops like 1/r, the outbound cloud 
will eventually reach a maximum velocity at the radius 
where these two forces are equal (which will occur at 
/3 = R C rit,c)- This allows us to compute a maximum 
outflow velocity for a given cloud: 

Vmax,c/Vcir — V^2 [(Rcrit.c 1) Rl(R’Crit.c)\ ^ (5) 

We overplot this equation on the data shown in Figure 
9 (using u m ax,c/2 to represent the predicted centroid of 
absorption-line tracing the outflow). The predicted and 
observed velocities are in satisfactory agreement for most 
of the data points. This is gratifying, given the simplicity 
of the model. In particular, the model allows us to under¬ 
stand the steep rise in v ou t/v C i r for relatively small values 
of R cr it,c (a regime where the ln(R cr it , c ) term is signif¬ 
icant), and a roll-over at large values where the model 
asymptotes to v ou t/v C i r — tj rna:r f 2u cv ; r = \Rcrit,d 2] // ■ 
We also note that Equation 4 predicts that clouds 
are rapidly accelerated. It is straightforward algebra to 
show that a cloud initially at rest at the starburst radius 
(/3 = 1) achieves half of v ma x,c defined by Equation 5 
at a radius /? < 4/3. This rapid radial rise in v ou t((3) is 
consistent with spatially-resolved maps of the outflows of 
ionized emission-line gas in nearby starbursts (e.g. Shop- 
bcll & Bland-Hawthorn 1998). 

4.1.2. Momentum-Driven Shell 

This model is similar to the one above, but assumes 
that the absorption arises in an expanding spherical shell 
of gas centered on the starburst. Just as in the case of 
the cloud, there is an outward force produced by the 


momentum supplied by the starburst, here given simply 
as F out = p*. There is also an inward force due to gravity, 
given by F in = M s v^ ir rJ 1 . We neglect pressure terms 
here. Defining f s to be the ratio of shell mass M s to 
the total mass interior to r, and then substituting for 
M s in the above equation we get F in = f s v^ ir /G. This 
expression for Fi n then defines the critical value for p* 
for the net force to be outward. In convenient units this 
is p C rit,s = 10 34 2 dynes for f s = 0.1 and v C i r = 100 km 
sec -1 . These values are listed in Table 2. 

In Figure 8, we plot this quantity vs. p* assuming that 
f s = 0.1. This plot is rather similar to the one for the 
cloud model, but here only 59% (23/39) of the galax¬ 
ies he above the critical value. We then plot v ou t /v C ir vs. 
P*/Pcrit,s in Figure 9. While this plot is qualitatively sim¬ 
ilar to that for the cloud, the scatter here is significantly 
larger (and the correlation is not statistically significant). 

Solving for the equation of motion for a shell we get: 

v s (/3)/v cir = V2[{p*/p C rit, s - l)ln(f3)\ 1/2 ( 6 ) 

In this case, the ratio of outward to inward force is inde¬ 
pendent of radius, and the shell does not reach a max¬ 
imum velocity. In Figure 9 we overplot the predicted 
value of v s /v C i r vs. the observed ratio v ou t/v C i r assuming 
that /3 = 2. The match to the data in this plot is sig¬ 
nificantly worse than for the cloud model. We conclude 
that this model is less successful than the cloud model 
in predicting both which galaxies will have outflows, and 
the observed outflow velocities. 

4.2. Further Considerations 

We concluded above that the momentum-driven clouds 
model does a better job in reproducing the properties of 
observed outflows, and now consider further aspects of 
this model. We found R cr it,c (the ratio of the momen¬ 
tum flux supplied by the starburst to the momentum 
flux necessary to overcome gravity) is a key parameter. 
We have then divided the sample into strong-outflows 
( Rcrit,c > 10) and weak-outflows ( R cr it,c < 10). 

To further understand these results, we return to our 
result in section 3 where we compared the momentum 
flux in the observed outflow to that injected by the star- 
burst. In Figure 7 we found a strong but non-linear cor¬ 
relation between p out and p*. The most forceful outflows 
(Pout > 10 34 dynes) usually carry a significant fraction of 
the starburst output (of order 100%). The less forceful 
outflows carry a much smaller fraction (of order 10%). 
These results can be understood in the following way. 
For an outflow that carries all the momentum flux sup¬ 
plied by the starburst we have M 0U tV 0 ut = 10 33 ’ 7 SFR 
dynes (where SFR is in M 0 per year). Dividing both 
sides by M out and converting SFR to cgs units, we get 
out = 760 SFR/M out km s -1 . The most forceful out¬ 
flows reach v out ~ 300 to 500 km s _1 , implying typical 
mass loading factors (M ou ,t/ SFR) of ~2 (consistent with 
the range shown in Figures 5 and 6). 

Figure 6 also shows that the mass-loading factor is ei¬ 
ther only weakly dependent on, or independent of, the 
bulk properties of the galaxy (M*, v C j r ), the starburst 
(SFR, r*), or the wind ( v ou t ). This would then imply 
that that the weak-outflows carry roughly the same (nor¬ 
malized) mass flux as the strong-outflows, they just do 
so at a significantly lower velocity and therefore carry 
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Fig. 8. — The rate of momentum supplied by the starburst plotted as function of the critical value of momentum flux needed to overcome 
gravity, and to hence drive an outflow (see text). The left panel (a) applies to a model of a population of clouds and the right panel (b) to 
a spherical shell of gas. The diagonal lines show p*/p C rit = 10, 1, and 0.1. For the cloud model roughly 95% of the galaxies lie above the 
critical value, compared to only 60% for the shell model. Outflows are observed in about 90% of the starbursts in our sample. The crosses 
represent the typical uncertainties (see Table notes). 
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Fig. 9. — The observed normalized outflow velocity ( v ou t/v c i r ) is plotted as a function of the ratio of the amount of momentum flux 
supplied by the starburst to the critical value needed to overcome gravity and drive an outflow (see text). The left (a) and right (b) 
panels show the respective cases for the cloud and shell models. The cloud model shows three regimes (indicated by dashed vertical lines). 
A ‘no-outflow’ regime (p* < p cr it ), a ‘weak-outflow’ regime in which v 0 ut/v c i r increases rapidly with increasing p*/p C rit, and a ’strong- 
outflow’ regime above p*/p cr it ~ 10 where v 0 ut/v c i r is typically 3 to 5. The over-plotted curve is derived from the equation-of-motion 
for a momentum-driven cloud, and agrees with the data in most cases. The right panel for the shell model is qualitatively similar but 
shows much greater scatter and does a poorer job of predicting which starbursts drive outflows. The plotted curve is derived from the 
equation-of-motion for the shell model, and is generally a poor fit to the data. See text. The crosses represent the typical uncertainties 
(see Table notes). 


a relatively small fraction of the momentum flux sup¬ 
plied by the starburst. This presumably reflects the fact 
that inward gravitational forces are significant in these 
cases. In Figure 10 we show that in nearly every case, the 
strong-outflows (ft* > 10p cr j t ) carry a significant fraction 
(of-order unity) of the momentum flux available from the 
starburst. This is not generally true for the other star- 
bursts. 

Finally, we can re-examine the strong empirical cor¬ 


relations of the observed normalized outflow velocity 
( Vout/vcir ) with both and SFR/area and SFR/M* (Fig¬ 
ure 2). So long as the ln(R cr it, c ) term in Equation 5 is 
insignificant (i.e. R cr it,c » 1), the expression can be 
rewritten as: 

( ft max ,c / U cir ) ^ 2 Rcrit,c — 2 P*/Pcrit,c ( 7 ) 

Since ft* oc SFR , ft C rit, c oc r*v^ ir , and v cir oc Mf 29 , 
equation 6 implies: 

Vmax,c/v c ir OC [(S F R/(r*Mft 58 )} 1 /2 (8) 
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Fig. 10.— The ratio of the momentum flux available from the 
starburst to the momentum flux in the observed outflow is plotted 
as function of the ratio of the starburst momentum flux to the 
critical value required to overcome gravity. This figure shows that 
most of the strong ’saturated’ outflows (p* > 10 p C rit - to the 
right of the dashed vertical line) carry a significant fraction of the 
available momentum flux, while the others generally do not. See 
text. The crosses represent the typical uncertainties (see Table 
notes). 

Given this expression, it is not surprising that the ob¬ 
served normalized outflow velocity correlates well with 
both SFR/rl and SFR/M * (Figure 3). 

4.3. Implications for ‘Sub-Grid Physics’ 

As described in section 1, both cosmological numerical 
simulations and semi-analytic models for galaxy evolu¬ 
tion often make use of simple prescriptions to deal with 
the effects of feedback from stars (‘sub-grid physics’). 
The current state of such models have been reviewed 
by Somerville & Dave (2015). We close the paper by 
briefly summarizing the implications of our results for 
these models. 

4.3.1. Comparing Data and Prescriptions 

We begin by comparing our data to predictions of out¬ 
flow velocities. One popular prescription is to assume 
that v out oc v C ir (e.g. Baugh et al. 2005; Bower et al. 
2006; Oppenheimer & Dave 2008; Dutton et al. 2010; 
Guo et al. 2011; Dave et al. 2013). It is clear from Fig¬ 
ures 1 and 2 that overall, this is a bad assumption: the 
ratio v ou t /u C j r spans about two orders-of-magnitude and 
has strong systematic dependences on both SFR/M* and 
SFR/area. Examination of Figure 9 shows that most of 
the scatter in the ratio of v ou t/v C i r is produced by the 
weak-outflows (R C rit,c < 10). With a few exceptions, 
the strong-outflows (R C rit,c > 10) have v ou t ~ 3 — 5 v C i r . 
We have explained above how these results can be un¬ 
derstood in the context of an outflow of a population of 
momentum-driven clouds. 

Next, we consider our results on the relationship be¬ 
tween mass-loading and outflow velocity. A common pre¬ 
scription is to assume that outflows use up the momen¬ 
tum flux supplied by the starburst, which results in a 


scaling such that SFR oc M 0 utv 0 ut■ This is equivalent 
to assuming that the mass-loading factor in the outflow 
M ou t/ SFR oc v~ut (e.g. Oppenheimer & Dave 2008; Dut¬ 
ton et al. 2010). Figure 11a compares this scaling rela¬ 
tion to the data, and shows it does not match the weak- 
outflows. As we have argued above, the scaling based 
on momentum conservation fails here because only the 
strong-outflows carry a substantial fraction of the avail¬ 
able momentum flux. The weak-outflows are slow be¬ 
cause they are ineffective at tapping this source of mo¬ 
mentum, not because they are strongly mass-loaded. The 
strong-outflows are largely consistent with this prescrip¬ 
tion. 

A related prescription is to assume the outflows carry 
a fixed fraction of the kinetic energy supplied by the star- 
burst, implying that the mass-loading factor scales like 
v~ut (e.g. Baugh et al. 2005; Somerville et al. 2008). Fig¬ 
ure 11a shows that the strong-outflows of typically carry 
about 10 to 50% of the available kinetic energy flux. This 
fraction can be much smaller in the weak-outflows in our 
sample. 5 The vast range in this fraction underscores 
the problem with adopting a simple proportionality be¬ 
tween the input and output kinetic energy fluxes for all 
outflows. 

Hybrid models with a mass-loading factor oc v~^ t at low 
velocities and oc v~ ut at higher velocities have also been 
adopted (Dave et al. 2013; Vogelsberger et al. 2013; 
Ford et al. 2015). These have disagreements with the 
data that are similar to those of the simpler momentum- 
based and energy-based prescriptions (severely failing for 
the weak-outflows). 

A different type of prescription assumes a fixed mass¬ 
loading factor ( M ou t oc SFR) and that the outflows all 
have the same velocity (e.g. Springel & Hernquist 2003; 
Schaye et al. 2010; Dave et al. 2013; Ford et al. 2015). 
As seen in Figure 11a, this prescription fails dramatically 
for the weak-outflows. We have shown that the outflow 
velocity depends on p*/p cr it in these cases. 

As noted above, most prescriptions assume that the 
outflow velocity is proportional to the circular velocity. 
Thus, the proportionality between the mass-loading fac¬ 
tor and the circular velocity in these models will be the 
same as that for the outflow velocity. This is shown more 
explicitly in Figure lib. While we observe a general in¬ 
verse correlation between the mass-loading term and cir¬ 
cular velocity (as is assumed in most prescriptions), it 
important to emphasize that the weak-outflows are gen¬ 
erally moving mass outward at velocities that are below 
the galaxy circular velocity (and certainly well below es¬ 
cape velocity). This is inconsistent with the assumptions 
in the model prescriptions. Considering only the strong- 
outflows (which are traveling at or above the escape ve¬ 
locity), we find a best-fit relation with a slope of -0.98 
and unit mass-loading factor at v C i r = 10 2 ' 25 km sec -1 . 
While this slope is consistent with a momentum-based 
scaling relation, the more direct test of such a relation 
is the plot of mass-loading vs. v out in Figure 11a. We 
also note that our empirical scaling relation is in good 

5 We note that this does not necessarily imply that most of 
the kinetic energy input from the starburst has been lost through 
radiative cooling. The bulk of this energy could reside in the fast 
(but tenuous) wind fluid whose ram pressure is driving the outflow 
seen in the absorption-lines (see section 1). 
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Fig. 11.— Left (a): A magnified version of Figure 6b showing the observed relation between the mass-loading factor in the outflow 
(M ou t/SFR vs. the outflow velocity v ou t). The plotted lines refer to various assumptions made in semi-analytic models or cosmological 
simulations (see text). The blue line shows the assumption that the outflow contains 100% of the momentum flux supplied by the starburst. 
The two red lines show the assumption that the outflow contains 100% (upper) and 10% (lower) of the kinetic energy flux supplied by 
the starburst. A hybrid ‘momentum-energy’ model is shown as a broken green line. The large cyan hollow circle assumes a constant 
outflow velocity and mass-loading factor. The blue data points are for the strong-outflows ( R C rit,c > 10) and the green points are for 
the weak-outflows ( R C rit,c < 10). Only for the strong-outflows do the assumed relations roughly match the data. Right (b): A magnified 
version of Figure 6c showing the observed relation between mass-loading and the galaxy circular velocity v c i r . The blue and red diagonal 
lines schematically indicate scalings often made in models and simulations (with mass-loading proportional to J; and to t?” J?, respectively. 
The blue and green symbols indicate the strong- and weak-outflows respectively. The dashed black line shows a least-square fit to the 
strong-outflows. See text for details. The crosses represent the typical uncertainties (see Table notes). 


agreement (for v C i r > 50 km sec -1 )with the recent FIRE 
simulations by Muratov et al. (2015) evaluated at z ~ 
0. These simulations explicitly calculate stellar feedback 
rather than adopting a prescription. 

In conclusion, some of the commonly adopted model 
prescriptions for outflow velocity and mass-loading agree 
with the observed properties of the strong-outflows. 
However, none of them adequately represent the prop¬ 
erties of our full sample. In particular, the assumption 
that the outflow velocities simply scale with the galaxy 
circular velocity is incorrect for the weak-outflows, and 
the assumptions about how the mass-loading scales with 
the outflow velocity lead to disagreements with the data 
that are over an order-of-magnitude in many weak out¬ 
flows. 

4.3.2. The Way Forward 

While current cosmological simulations with increas¬ 
ing numerical resolution are improving the capability to 
model outflows in a more realistic fashion (e.g. Cen 2014; 
Vogelsberger et al. 2013; Schaye et al. 2015; Henriques 
et al. 2015; Mitra et al. 2015), high-resolution simula¬ 
tions of individual galaxies that are able to calculate ab 
initio more of the relevant physical processes are clearly 
needed (e.g. Cooper et al. 2008, 2009; Shen et al. 2013; 
Hopkins et al. 2013; Muratov et al. 2015). Ideally, these 
simulations should be closely compared to (and influ¬ 
enced by) both the results in this paper and to detailed 
multi-phase observations of actual galactic winds in the 
local universe. 

In the meantime, a formulation based on the results 
in this paper could be useful. That is (given a star- 
formation rate, a radius of the star-forming region, and 


a halo circular velocity) calculate p*/p C rit,c ( Rcrit,c )• Use 
this to assign an outflow speed using our equation 5 above 
(and see Figure 9). It may be case that most stars over 
cosmic history form in systems with R cr it,c > 10 (lead¬ 
ing to strong-outflows). In this case, the simple scal¬ 
ing in Equation 7 could be used: v ou t = v maXiC /2 ~ 

Vcir \j Rcrit.c/2“ 

The mass-loading term as defined in mod¬ 
els/simulations is only really applicable to the strong- 
outflows, which are the only cases in which the outflow 
velocity typically exceeds the galaxy escape velocity. 
For these, the best-fit scaling relation for mass-loading 
shown in Figure lib could be adopted. However, we 
emphasize that actual values of the mass-loading term 
were derived from our data under the assumption that 
the (absorption-weighted) size of the outflow is twice 
the radius of the starburst. A different scale factor will 
result in a different constant of proportionality between 
mass-loading and circular velocity for strong-outflows. 

5. CONCLUSIONS 

We have analyzed outflows as traced by far-UV ab¬ 
sorption lines for a combined sample of 39 low-redshift 
starburst galaxies. The majority of these starbursts (24) 
have properties that make them excellent local analogs 
to typical high-redshift Lyman Break Galaxies. This 
is in contrast with most previous investigations of low- 
z outflows, which have analyzed far-IR-selected star- 
bursts using the optical Na I absorption-line. Com¬ 
pared to most studies of higher-redshift UV-selected 
galaxies, our individual spectra are higher in S/N and 
the sample spans much wider ranges in star-formation 
rates (SFR), galaxy stellar masses (M*), circular ve- 
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locities (vcjr), metallicities, specific star-formation rates 
(SFR/M*), and SFR/area. 

We begin by summarizing the most robust empirical 
conclusions. These all involve directly measured quan¬ 
tities without uncertain assumptions. We found that 
the outflow velocity of the ionized gas seen in absorp¬ 
tion correlates only weakly with the galaxy stellar mass 
(M*) and the galaxy circular velocity (rw), but shows 
significant correlations with SFR/M*, and (especially) 
with SFR and SFR/area. Notably, the outflow velocity 
in our sample increases with increasing SFR/area, but 
appears to ‘saturate’ at roughly 300 to 500 km s _1 for 
SFR/area > 6M 0 kpc -2 . The ratio of v ou t/v C i r spans 
nearly two orders-of-magnitude and correlates strongly 
and positively with both SFR/M* and SFR/area. This 
ratio reaches typical values of ~ 3 to 5 for starbursts with 
high SFR/M* and SFR/area, well in excess of the galaxy 
escape velocity. 

Our conclusions regarding outflow rates are more un¬ 
certain. These require values for the total gas column 
density and characteristic radius of the outflow as seen 
in absorption. We have used weak (optically-thin) metal 
absorption-lines to find that the column density of ion¬ 
ized gas is over an order-of-magnitude larger than of the 
neutral gas, and to show that the total implied hydro¬ 
gen column density of-order 10 21 cm~ 2 . We simply pa¬ 
rameterized the effective (absorption-weighted) radius of 
the absorbing gas as a multiple of the starburst radius 
[r 0 ut = /3r*). With a choice /3 = 2, the typical esti¬ 
mated mass outflow rates range from ~ 1 to 4 times the 
SFR. The ‘mass-loading’ factor (M out /SFR) correlates 
only weakly and inversely with SFR, v C ir , or M*, and 
shows no correlation with v 0 ut or r*. 

We then considered two simple analytic models for the 
outflow: a population of clouds accelerated by the net 
sum of gravity and the momentum-flux supplied by the 
starburst, and a similar model in which the clouds were 
replaced by a spherical shell of gas. We compared the 
momentum flux provided by the starburst (p*) to the 
critical value needed to overcome gravity ( p cr it ) (which 
is different for the cloud and shell models). On this basis, 
the cloud model did a significantly better job of predict¬ 
ing both which galaxies should have outflows and the 
outflow velocities themselves. 

Further consideration of the model of momentum- 
driven clouds led to some interesting insights. We found 
that the ratio ( R cr it,c = P*/ Pcrit,c) determines the ob¬ 
served outflow velocity. We identified three regimes: a 
‘no-outflow’ regime where R cr it,c < 1, a ‘weak-outflow’ 
regime in which v 0 ut/v C ir increased with increasing 
Rcrit,a and then a ‘strong-outflow’ regime ( R cr it,c > 10) 
in which Vout/v C ir ~ 3 to 5. We showed that this be¬ 
havior can be rather well-represented by the equation- 
of-motion derived for this simple model of momentum- 
driven clouds. The strong-outflow regime corresponded 
to cases in which the estimated momentum flux in the 
outflowing gas ( p ou t ) is comparable to the total momen¬ 
tum flux supplied by the starburst (/>*), and in which the 
outflows are faster than the galaxy escape velocity. In the 
weak-outflow and no-outflow cases it is likely that grav¬ 
ity plays a significant role in the cloud dynamics and the 
coupling of the starburst momentum flux to outflowing 
gas is ineffective. 


We showed that none of the simple prescriptions used 
to capture the essential features of winds in typical cos¬ 
mological simulations and models of galaxy evolution 
agree with the empirical results for the full sample we 
have studied. In particular, there is no tight, simple scal¬ 
ing observed between v ou t and v C i r (the ratio of the ve¬ 
locities spans nearly two orders-of-magnitude), and v ou t 
varies by over a factor of 50 (inconsistent with the so- 
called ‘constant wind’ model). We see no evidence that 
the mass-loading factor (M out /SFR) correlates inversely 
with Vout, as is sometimes assumed. All these disagree¬ 
ments are severe for the weak-outflows. For most of the 
strong-outflows, at least some of the prescriptions can 
successfully describe the data within the uncertainties in 
our derived parameters. 

We believe the empirical results in this paper can be 
understood on the basis of a simple physically-motivated 
model of a momentum-driven outflow of a population of 
clouds. Our results could be used to define a simple pre¬ 
scription for the ratio of outflow to circular velocity based 
on the value of R cr it,c (which can be estimated from the 
model/simulated galaxy star-formation rate, the radius 
of the star-forming region, and the galaxy circular veloc¬ 
ity). The mass-loading term can also be parameterized 
as being proportional to v~ ir , but only for the strong- 
outflows ( R C rit,c > 10). 

In closing, we emphasize that our hope is that the re¬ 
sults presented in this paper will provide helpful guidance 
as to the treatment of galactic winds in future numerical 
simulations and theoretical models, as well as providing 
new insights into the nature of galactic winds driven by 
massive stars. 
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TABLE 1 

Description of Galaxy Properties. 


Galaxy 

z 

(kpc) 

Log M* 
(Log M 0 ) 

V c . 

cir 

(km s _1 ) 

SFR d 

(Moyr- 1 ) 

Log SFR/M* 
(log yr- 1 ) 

Log SFR/area 

(M 0 y r_1 kpc -2 ) 

12 + log[0/H] e 

T f 

T uv 

Sample 

J002101+005248 

0.09840 

0.53 

9.3 

83 

15 

-8.1 

0.93 

8.20 

0.8 

cos 

J005527-002148 

0.16744 

0.32 

9.7 

108 

24 

-8.3 

1.56 

8.28 

2.1 

cos 
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0.14668 

1.37 

10.3 

161 

37 

-8.7 

0.50 

8.39 

2.8 

cos 

J021348+125951 

0.21902 

0.39 

10.5 

184 

19 

-9.2 

1.30 

8.74 

1.9 

cos 

J080844+394852 

0.09123 

0.08 

9.8 

115 

8 

-8.9 

2.33 

8.74 

1.2 

cos 

J082354+280621 

0.04722 

0.34 

8.6 

52 

10 

-7.6 

1.12 

8.23 

2.4 

cos 

J092159+450912 

0.23499 

0.78 

10.8 

225 

29 

-9.3 

0.89 

8.67 

1.9 

cos 
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0.18072 

0.69 

9.1 

72 

10 

-8.1 

0.54 

8.09 

0.4 

cos 
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0.67 

9.4 

88 

11 

-8.4 

0.60 

8.19 
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cos 
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77 
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-8.3 
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cos 
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0.33 
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-8.7 
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cos 
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9.6 

94 
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cos 
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9 
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48 
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0.31 
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0.19 
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23 
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NGC 1705 

0.00211 

0.24 

8.6 

55 

0.16 

-9.4 

-0.35 

8.07 

0.15 

FUSE 
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172 

21 

-9.1 

0.44 

8.40 

0.9 

FUSE(LBA) 

M 83 

0.00172 

3.20 

10.7 

210 

3.5 

-10.2 

-1.26 

8.85 

2.6 

FUSE 

NGC 5253 

0.00135 

0.18 

9.3 

83 

0.33 

-9.8 

0.21 

8.15 

1.2 

FUSE 

IRAS 19245-4140 

0.00945 

0.42 

9.1 

72 

2.1 

-8.8 

0.29 

8.22 

0.45 

FUSE 

NGC 7673 

0.01131 

1.80 

10.0 

132 

4.8 

-9.3 

-0.63 

8.45 

1.4 

FUSE 

NGC 7714 

0.00935 

0.50 

10.2 

151 

6.9 

-9.4 

0.64 

8.47 

3.6 

FUSE 


The half-light radius of the starburst as measured with UV imaging (see text for details). The measurement uncertainties are typically <0.1 dex. 

The total galaxy stellar mass derived from near-IR and/or multiband optical photometry. Based on comparisons to the measured dynamical masses (Overzier et al. 2011), the typical uncertainties 
are ±0.3 dex. 

C The characteristic circular velocity based on the fit between M* and v c i r in Simons et al. (2015). The scatter in this relationship is ± 0.1 dex. 

The star-formation rate based on UV and IR photometry (see A15). A standard Kroupa/Chabrier IMF is assumed. Comparisons to SFRs derived from the extinction-corrected Ha luminosity 
imply an uncertainty of ± 0.2 dex. 

e The oxygen abundance (12 + log[0 / H]) derived from the Pettini Sz Pagel (2004) calibration of the 03N2 strong-line method. These authors quote an uncertainty of ±0.14 dex. 

f 

The effective optical depth in the FUV defined as —In^Luy^bs / Luv,int (see text). The uncertainties are caused by systematic errors in Luv,int as derived from UV plus IR photometry. We 
estimate the resulting uncertainty on tuv to be ~ ± 0.2 


Systematic Properties of Galactic Winds 






TABLE 2 

Description of Outflow Properties. 


Galaxy 

Kut 

(kms -1 ) 

Mjut 

(M 0 yr x ) 

Log M 0 ut/SFR 

Log PSut 
(log dynes) 

Log 

(log dynes) 

Log P crit,c 
(log dynes) 

Log l4it, s 

(log dynes) 

Log N h 
(log cm -2 ) 

J002101+005248 

320 

33 

0.35 

34.8 

34.9 

33.4 

33.9 

21.0 

J005527-002148 

410 

26 

0.04 

34.8 

35.1 

33.4 

34.3 

21.3 

J015028+130858 

360 

97 

0.41 

35.3 

35.3 

34.4 

35.0 

21.3 

J021348+125951 

500 

39 

0.31 

35.1 

35.0 

33.9 

35.3 

20.8 

J080844+394852 

560 

9.0 

0.02 

34.5 

34.6 

32.8 

34.4 

20.6 

J082354+280621 

510 

34 

0.55 

35.0 

34.7 

32.8 

33.1 

21.4 

J092159+450912 

480 

74 

0.40 

35.4 

35.2 

34.4 

35.6 

20.9 

J092600+442736 

350 

48 

0.66 

35.0 

34.7 

33.4 

33.6 

20.8 

J093813+542825 

280 

37 

0.52 

34.8 

34.7 

33.5 

34.0 

21.0 

,1102548+362258 

240 

30 

0.60 

34.7 

34.6 

33.4 

33.7 

20.7 

J111244+550347 

450 

30 

0.02 

34.9 

35.1 

33.7 

34.9 

21.0 

J111323+293039 

450 

99 

1.14 

35.4 

34.5 

33.8 

34.1 

21.0 

,1114422+401221 

290 

45 

0.70 

34.9 

34.6 

33.9 

34.6 

20.9 

J141454+054047 

30 

3.5 

-0.16 

32.8 

34.4 

33.0 

32.9 

21.0 

J141612+122340 

400 

15 

-0.19 

34.6 

35.0 

33.3 

34.7 

21.0 

,1142856+165339 

330 

47 

0.53 

35.0 

34.8 

33.6 

34.1 

20.8 

J142947+064334 

360 

21 

-0.11 

34.6 

35.1 

33.2 

34.7 

21.3 

J152141+075921 

280 

21 

0.55 

34.6 

34.5 

33.3 

34.1 

20.4 

,1152521+075720 

350 

35 

0.59 

34.9 

34.6 

33.4 

34.0 

21.0 

,1161245+081701 

450 

28 

-0.11 

34.9 

35.2 

33.5 

34.7 

21.0 

J210358-072802 

500 

46 

0.05 

35.2 

35.3 

34.2 

35.7 

21. 0 

Haro 11 

160 

45 

0.10 

34.7 

35.2 

34.3 

34.9 

21.3 

VV 114 

130 

60 

-0.04 

34.7 

35.5 

34.9 

35.6 

21.1 

NGC 1140 

34 

4.8 

0.76 

33.0 

33.6 

33.5 

34.0 

20.9 

SBS 0335-052 

38 

2.3 

0.86 

32.7 

33.2 

32.2 

32.1 

21.0 

Tol 0440-381 

120 

33 

0.82 

34.4 

34.4 

34.2 

34.7 

20.9 

NGC 1705 

21 

1.0 

0.80 

32.1 

32.9 

32.7 

33.2 

20.4 

NGC 1741 

38 

4.6 

-0.12 

33.0 

34.5 

33.6 

34.3 

21.1 

I Zw 18 

< 10 

< 1 

< 1.80 

< 31.8 

31.9 

32.1 

31.3 

21.0 

NGC 3310 

180 

22 

0.90 

34.4 

34.1 

33.7 

34.4 

21.1 

Haro 3 

< 10 

< 1.0 

< 0.39 

< 31.8 

33.3 

33.2 

33.5 

21.1 

NGC 3690 

100 

12 

-0.52 

33.9 

35.3 

34.3 

35.7 

21.3 

NGC 4214 

34 

1.0 

0.89 

32.3 

32.8 

32.6 

33.5 

21.1 

Mrk 54 

58 

13 

-0.21 

33.7 

35.0 

34.3 

35.1 

20.8 

M 83 

< 10 

< 6.4 

< 0.26 

< 32.6 

34.2 

34.9 

35.5 

20.9 

NGC 5253 

17 

0.6 

0.26 

31.8 

33.2 

33.2 

32.9 

21. 2 

IRAS 19245-4140 

65 

5.4 

0.41 

33.3 

34.0 

33.1 

33.6 

20.7 

NGC 7673 

82 

30 

0.80 

34.2 

34.4 

34.3 

34.7 

21.0 

NGC 7714 

46 

4.6 

-0.18 

33.1 

34.5 

33.9 

34.9 

21.4 


a The outflow velocity of the strong interstellar absorption-lines arising in the warm ionized gas. The uncertainties are generally dominated by systematic uncertainties in the fit to the adjacent 
continuum. A15 estimate the resulting uncertainties to be ± 0.05 dex. 

k The mass outflow rate derived from the outflow velocities, an estimate of the outflow column density, and the column-density-weighted mean radius of the absorbing gas (see text). A uniform 
column density of 10 21 cm -2 was assumed, leading to an estimated uncertainty of ± 0.25 dex (see Figure 4). A uniform radius of 2 r* was assumed for the calculation. 

C The momentum flux in the outflow. This is just the mass outflow rate (Column 3) times the outflow velocity. The uncertainties will be dominated by the former. 

The momentum flux from the starburst. This is assumed to be proportional to the star-formation rate (see text). This conversion may introduce systematic uncertainties that are hard to quantify 
(Leitherer et al. 1999). 

e The critical momentum flux needed from the starburst to overcome gravity acting on a cloud with a column density of 10 21 cm -2 located at the radius of the starburst. The uncertainties are 
all systematic, driven by these adopted values. 

f 

The effective column density of dusty gas along the line-of-sight, based on the effective FUV optical depth (Column 10, Table 1), the Calzetti et al. (2000) FUV starburst attenuation law, and 
the assumption of a constant dust/gas-phase-metals ratio. The uncertainties are similar to the observed spread in the distribution (± 0.3 dex). 
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